Map of metastable states for thin circular magnetic nano-cylinders 
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Nano-magnetic systems of artificially shaped ferromagnetic islands, recently became a popular 
subject due to their current and potential applications in spintronics,— magneto-photonics^ and 
superconductivity k When the island size is close to the exchange length of magnetic material (around 
15 nm), its magnetic structure becomes markedly different.—^ It determines both static and dynamic 
magnetic properties of elements, but strongly depends on their shape and size.— ^ Here we map this 
dependence for circular cylindrical islands of a few exchange lengths in size. We outline the region of 
metastability of "C"-type magnetic states, proving that they are indeed genuine and not a result of 
pinning on particle imperfections. A way to create the smallest particles with guaranteed magnetic 
vortex state at zero field becomes evident. It is expected that the map will help focus the efforts in 
planning of experiments and devices. 
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Since the first theoretical model for magnetic vortex 
in circular cylinder, taking into account long-range mag- 
netic dipolar interaction^ and experimental discovery of 
magnetic vortex core^ there have been a number of devel- 
opments in this area. The question of magnetic ground 
states in cylinders of various sizes was settled^ 6 ^ and 
thoroughly tested 7 - both experimentally and numerically. 

Yet, the most important property of ferromagnets, re- 
sponsible for wealth of their present and future appli- 
cations, is the ability to be in one of many metastable 
states. Fortunately, the number of distinct metastable 
states for nano-magnets is much smaller than in ordinary 
ones, making it feasible to enumerate and map them all. 
Such a map for circular cylinders in coordinates of R/Le 
vs L/Le, where L and R are the particle's thickness and 
radius, respectively, and Le = \/ AirC / /io7B-<Wj ls the 
exchange length, is shown in Figure. C and Ms are the 
exchange constant and the saturation magnetization of 
material, respectively, fiQ is the permeability of vacuum 
in SI units or 1 in CGS and 73 is 1 in SI or 4-7T in CGS. 
Let us discuss the map first, before diving into the details 
on how the additional lines were derived. 

The simplest magnetic states of nano-magnets 
are (quasi-)uniform ones: in-plane and out-of-plane. 
Neglecting quasi-uniformity, making only a small 
correction^ the energies of these states^ are equal in 
cylinders with aspect ratio g = L/R = gj rj 1.8129523. 
At this critical aspect ratio (straight line "a" in Figure) 
there is no energy barrier between the states, and the 
equilibrium line coincides with the stability line for each. 

Solving for geometries with the energy of the vortex 
state equal to the energies of the (quasi-)uniform states 
(in-plane or out-of-plane, depending on aspect ratio), 
yields^ line "b" . Lines "a" and "b" separate the regions 
with different ground states of the particle, which can be 
reached experimentally, e.g., by annealing or thermal re- 
laxation simply through waiting for a sufficient time. 

Study of the vortex stability with respect to the dis- 
placement from the center of circular face produced 6 line 



"c" on the map. Below it (at a low R/Le) the vortex 
state is impossible without external stabilization. This 
is the lower bound for the region of coexistence of the 
vortex and the (quasi-) uniform states. An upper bound 
of another metastability region, where the vortex state 
coexists with the out-of-plane uniform oneii is shown as 
line "d," above which the uniform state is unstable. From 
general considerations the line "b" should coincide or lie 
below the line "d" . It is not so (by hardly experimentally 
detectable amount) due to neglected magnetization vari- 
ation along the cylinder thickness, as discussed in RefflTI. 

Here we fully outline the metastability regions by in- 
vestigating the bending instability of the quasi-uniform 
"leaf" state (line "e"), as well as stability of vortex with 
respect to the core expansion (line "g"). It was puzzling 
at first to obtain the stability line "e," crossing the equi- 
librium line "b," but it turned out that destabilization 
of the "leaf" state does not imply the vortex formation. 
Another energy minimum appears, capturing the magne- 
tization in a half-bent "C" state. Studying its stability 
yields line "f" , where the capturing minimum vanishes 
and only the vortex state remains stable. 

In-plane magnetization of "C" and "leaf" is unstable 
in elements with g < gj, giving rise to a "dip" of the 
vortex-only region down to the tricritical point of the 
diagram. If the in-plane state is created there (e.g., by 
a transient application of magnetic field), it turns into 
an out-of-plane state (because g < gj) and collapses into 
vortex (since the region is above the line "d" ) . 

A number of experiments and numerical result s 7 ' 12 ! 13 
confirm this map. Let us present only few, supporting the 
new lines most directly. First, a "unfitting"— point from 
an experiment fi 2 - where a remanent in-plane magnetiza- 
tion state was found in dots well above the equilibrium 
line "b." It fits precisely in the metastability region for 
the "leaf" state. Even more direct confirmation comes 
from in-plane hysteresis measurements,— where in-plane 
hysteresis loops of supermalloy dots of different sizes were 
classified. The data show sketchily the boundary, sepa- 
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FIG. 1: Map of metastable circular nano-cylinders. Solid lines represent the equality of the equilibrium energies of the bordering 
ground states. Dashed lines correspond to the stability loss. More than one state (represented by symbols with arrows) can 
be stable in the shaded regions. The leftmost in each group of symbols indicates the ground state. The formula gives the 
asymptotic to the "C"-vortex transition line, and the other dash-dot line is the equilibrium radius of the vortex core. 



rating the smaller dots with narrow rotational hystere- 
sis, and bigger dots, where the loop is of nucleation type. 
It is exactly the stability line "f" . Below the line, re- 
duction of the field from in-plane saturation creates the 
"C" magnetization state, which is then rotated (as field 
passes zero) and straightened to become the "leaf" again. 
Above the line, the "C" state is unstable and vortex is 
nucleated instead, which must be annihilated, produc- 
ing nucleation-type hysteresis loops. All the boundary 
points^ fit the map, but only three of them, most tightly 
converged around the line "f" are plotted. 

Let us now briefly describe the steps to obtain the lines. 
The starting point is provided by a recently developed 
approach^ to systematic generation of Ritz functions for 
magnetization distributions in thin nano-cylinders and 
its generalizationi^ii^ We will not follow it here in detail, 
but only present the resulting function as an anzats 

where z — X/R + iY/R is the normalized complex co- 
ordinate on the face of the circular cylinder, c is real, 
\z\ < 1, and a, p are (in general) complex Ritz parame- 
ters, whose values are determined by minimizing the total 



energy of the particle (consisting of the exchange and the 
magnetostatic contributions). Normalized magnetization 
vector m — M /Ms (\m\ = 1) components are expressed 
as mx + vm,Y — 2w/(l + ww), mz = (1 — ww) /(I + ww), 
where a complex (but not holomorphic) function w(z,z) 
is equal to f(z) if |/(z)| < 1 and f(z)/\f(z)\ if f(z) > 1. 

All the variety of magnetization distributions in nano- 
scale circular dots is given by ([1]) with different values of 
c, p and a. When p = 1 , the magnetization has no normal 
components to the particle's circular side. If \a\ < c, it 
is a displaced magnetic vortex. When a — 0, the vortex 
is centered (the case studied by Usov and Peschanyi to 
find the vortex core radius Ry = R/c). Consideration of 
< \a\ -C 1 yields^ line "c" in Figure. When c = and 
a = 0, the magnetization is uniform out-of-plane (mz — 
1, and mx = my = 0). Expanding its energy at < 
c«l allows to find the cylinder radius at which the out- 
of-plane configuration loses stability and is transformed 
into vortex.— The case of \p\ > 1, c = and \a\ > 
corresponds to the quasi-uniform "leaf" configuration. 9 
When p is imaginary, the configuration is changed into 
"flower," which is not realized with isolated dots, but can 
be stabilized by interaction in dot arrays*^ 

In finding the lines "e" and "f" we permit non-zero 
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values of c in the quasi- uniform state (p > 1 and \a\ > c), 
bending the "leaf" . The equation |T|) is parametrized by 
letting c = xa, where \x\ < 1, assuming that p > 1, a 
is real and a > ao, such that \f(z)\ > 1 always. The 
magnetization is independent of a and at x = coincides 
with that in Ref M (noting that x of Ref M is called p here) . 
Evaluation of magnetostatic energy of this distribution is 
cumbersome (we use the magnetic charges formalism and 
factorization)^ Normalized to ^q^bM^LttR 2 up to the 
fourth order in 1/p, it is 
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The functions e", e 2 , and e 4 are exactly the same as in 
Refi (to which expression this new result passes in a 
limit), and the three new functions are 

oc 

eM= f^m^m dk<0 , 



, , jFjkg) k 2 (J 2 (k)- J|(fc)) + 8J 1 (fc)J 3 (fc) „ ^ n 
e 42 (g) = / — dfc > 0, 



/F(fcg) Jj{]f) + 8J 2 (fc) - 9J|(fc) 
e44(5) = - / — £ g dfc < 0, 



where F(a;) = 1 — [1 — cxp(— x)]/x and J,(x) are the 
Bessel's functions of the first kind. They can be expressed 
in closed form through the complete elliptic integrals. 

The exchange energy, calculated as (C/2) ^(Vmi) 2 , 
in the same normalization and order in 1/p is 
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where p — R/Le- 

The total energy e 1 = e EX + e\fg is a polynomial in 
1/p and It can be easily minimized analytically. The 



solution for equilibrium p always exists and stable. For 
x, surprisingly, the equations gave two solutions. In cer- 
tain geometries the solution x — (the "leaf" state) is 
stable, whereas in other geometries the energy minimum 
at x — turns into the maximum, and another solution 
at \x\ < 1 (which is a "C"-type state, or a bent "leaf") 
becomes stable. For even larger cylinders the only stable 
solution has x > 1, which is vortex. The critical lines "e" 
and "f" result from simultaneously solving three equa- 
tions de/dp = 0, de/dx — and d 2 e/dx 2 — 0. These 
equations are polynomial and the critical lines can be 
expressed in a closed form. 

Consideration of vortex core expansion instability re- 
quires to allow for the core to be bigger than the par- 
ticle size (the case not considered in Refs.4.8). That is 
p = 1, a = and c < 1 (or c = cos a, a £ [0, ir/2} ) 
in Eq. ([T]). Only the surface magnetic charges, m z (r) — 
2/(1 + c 2 r 2 ) — 1, are present. Their normalized magne- 
tostatic energy (the self-energy of both faces and their 
interaction) can be written as e 2 MS = gu(0, ot) — gu(g, oc). 
The function gu is factored as 



gu(h,a) = I < 
o 
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The exchange energy after direct integration is 
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The critical line "g" was calculated numerically, by solv- 
ing a system of equations: de 2 / da = and d 2 e 2 / da 2 = 
0, where e 2 = e% x + e 2 MS . 
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